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ABSTRACT 

Equations describing the temporal and spatial behavior of the kinematic moisture and heat flux are introduced 
in this study. In these nonlinear equations, the contribution by diabatic processes to the large-scale flux is 
composed of two parts. One part is associated with a Rayleigh damping term, while the other arises from 
temporal and spatial changes in the pressure gradient term. 

The influence of diabatic processes on the large-scale moisture fluxes depends greatly on the degree of balance 
between forcing and damping terms in the governing equations. The existence of a near balance requires a 
reduction in the large-scale horizontal geostrophic wind speed. From a scale analysis of the moisture flux equations 
it is argued that reductions in the large-scale horizontal wind speed, observed within major cumulus cloud 
systems, help conserve large-scale moisture fluxes. The deviation of the wind from geostrophic conditions is 
easily estimated. This wind modification induces secondary vertical circulations that contribute to the convergence, 
creating or supporting long-lived mesoscale flows. In the tropics the wind modification has an antitnptic rela- 

These diagnostic findings suggest possible modifications to the wind field in the application of cumulus 
parameterization, and may be important in diabatic initialization of numerical weather prediction models. 


1. Introduction 

In the extratropics, the large-scale atmospheric flow 
fluctuates about known balance conditions. Away from 
the major centers of cyclonic and anticyclonic activity, 
the horizontal wind components closely approximate 
geostrophy in the free atmosphere. Within strong cy- 
clonic centers the flow is modified by curvature and 
by diabatic effects. To forecast the wind field in the 
presence of diabatic activity would require solutions 
of the complete equations of motion with physical pro- 
cesses included, which is not always practical. Other- 
wise, to get a snapshot of the instantaneous horizontal 
wind field requires solving the nonlinear diabatic bal- 
ance equations. The simple wind laws, like the gradient- 
wind approximation, do include the effect of curvature 
and some indirect consequences of diabatic heating. 
This study outlines a simplified approach derived from 
first principles that identifies how diabatically induced 
modifications in the horizontal wind components in- 
fluence the moisture and sensible heat fluxes. This 
study provides necessary additional information about 
the effect of diabatic processes on the atmospheric dy- 
namics for use in model initialization and cumulus 
parameterization. 

Many years of observations have shown that diabatic 
processes can significantly enhance cyclonic develop- 
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ment and influence the general circulation. Many 
journal articles and textbooks testify to this. A small 
subset includes Petterssen ( 1956), Riehl and Malkus 
( 1958), Kuo ( 1965, 1974), Yanai et al. ( 1973), Riehl 
(1979), Donner et al. (1982), Houze and Hobbs 
(1982), and Konig and Ruprecht (1989). Modeling 
studies also have verified that moist processes are im- 
portant for rapid development (Kuo and Reed 1988; 
Reed et al. 1988). Additionally, diabatic initialization 
techniques are used to enhance cyclonic circulations 
and low-level convergence to help speed up precipi- 
tation processes in numerical forecast models (e.g., 
Tarbell et al. 1981; Krishnamurti et al. 1983; Krish- 
namurti et al. 1984; Donner 1988; Wang and Warner 
1988; Turpeinen et al. 1 990; Daley 1991). 

There also are many studies on the effect of clouds 
on the large-scale dynamics, including Malkus ( 1952), 
Schneider and Lindzen ( 1976), Sui and Yanai (1986), 
Cho (1985), and Konig and Ruprecht (1989). The 
vertical transport of atmospheric properties in clouds 
is important because the cloud vertical velocity is much 
larger than that found in the environment, especially 
in the tropics. In contrast, the magnitudes of the hor- 
izontal flow in the cloud and environment are of the 
same order, but according to the above studies, mixing, 
vertical transport, and other processes are responsible 
for some deceleration of the mean flow by “cumulus 
friction” (Malkus 1952; Schneider and Lindzen 1976). 
Several articles propose parameterizations of cumulus 
convection effects on the large-scale vorticity or mo- 
mentum. Approaches include mass entrainment and 
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detrainment concepts (Ooyama 1971 ), vertical trans- 
port and properties from simple cloud models (Shapiro 
and Stevens 1980), the introduction of drag terms (Es- 
bensen et al. 1987; Austin and Houze 1973), and other 
representations of cumulus friction (Schneider and 
Lindzen 1976). In this study a different method is pro- 
posed to estimate the cumulus friction process. 

Gill ( 1980) has examined tropical circulation forced 
by localized heating, while Heckley and Gill (1984) 
have studied transient features in the tropics when the 
heating is impulsively started. In these studies the role 
of diabatic heating on the large-scale flow is modeled 
by adding a forcing term to the continuity equation in 
the linearized shallow-water model. The forcing mod- 
ifies the pressure and its gradients, which in turn modify 
the horizontal wind field. This model is conceptually 
simple yet captures the essence of the large-scale re- 
sponse to diabatic forcing arising from thermal changes. 
Gill (1980) heuristically added a Rayleigh damping 
term with constant coefficient over the entire domain 
of his linear model, which produced a steady state, but 
this was not intended to represent the cumulus friction 
process, and it provided no insight on why cumulus 
friction occurs, its region of influence, or consequences. 
Holton and Colton (1972) in their linear vorticity 
model also found that damping was needed to obtain 
agreement with observations. 

Lorenz (1978) has identified some optimal atmo- 
spheric energetic arrangements, but local optimal re- 
lationships remain unknown if they exist at all. One 
can speculate that cumulus friction is partially an at- 
mospheric response for conserving basic large-scale at- 
mospheric properties or optimizing fundamental pa- 
rameters like the kinetic energy, enthalpy, etc. Clearly, 
temporal and advective changes in moisture flux, heat 
flux, and entropy are influenced by cumulus friction. 

In this diagnostic study it is found that mesoscale 
moist processes introduce into the moisture and heat 
flux equations a Rayleigh damping term that depends 
upon the condensation and evaporation rate. This fric- 
tion allows changes in the large-scale moisture fluxes 
to remain relatively small provided a near balance exists 
between the forcing and damping terms. Scale analysis 
shows that this balance involves the pressure gradient, 
Coriolis, and Rayleigh damping term. Balance further 
requires a reduction in the large-scale horizontal geo- 
strophic wind speed, that is, the creation of an ageo- 
strophic wind. Clearly, if the horizontal wind is reduced 
excessively, then the balance condition breaks down. 
Our hypothesis is that for long-lived mesoscale con- 
vection the large-scale horizontal winds do not differ 
greatly from that needed to establish a near balance 
between the pressure gradient, Coriolis, and Rayleigh 
damping terms. It is found diagnostically and in nu- 
merical applications that the ageostrophic wind en- 
hances the existing convergence and establishes a me- 
soscale vertical circulation within the large-scale flow 
to help balance the modified horizontal wind. 


To illustrate our hypothesis we examine the behavior 
of the kinematic moisture and heat flux. The derivation 
of these equations is quite simple, but they yield sur- 
prising information. In these equations the contribu- 
tion of diabatic heating to the horizontal wind com- 
ponents is naturally partitioned into two parts. One 
diabatic contribution is contained within the geo- 
strophic wind component (pressure gradient term), 
while the remaining one is associated with the Rayleigh 
damping (condensation ) or forcing (evaporation ) term. 
To examine the importance of the terms in the flux 
component equations, a scale analysis is preformed and 
simplified solutions are found. For nontropical regions 
our approximate solutions are given in terms of the 
geostrophic flow. In some of our illustrative examples 
a time scale less than the inverse of the Coriolis param- 
eter is required; however, the flux equations are valid 
even when this is not true. In tropical regions any ten- 
dency toward moisture flux conservation enhances the 
relationship between the pressure gradient term and 
the Rayleigh damping term. This induces secondary 
antitriptic circulations. A simple hurricane model is 
presented to illustrate the wind modifications associated 
with any enhancement toward moisture flux conser- 
vation. Also tested are applications in cumulus param- 
eterization. 


2. The equations 


The inviscid equations (Haltiner and Williams 1980, 
pp. 16, 17, and 308 ) for large-scale flow describing hor- 
izontal wind components u and v y the temperature T y 
potential temperature 6 y the mixing ratio q, and pres- 
sure p are 


du 

~di 

dv 
dt : 

dd 
dt : 

dq 

dt 


1 dp r 
p dx 

0) 

1 dp 

-T fu, 

P dy 

(2) 

L6Q 
Tc P ’ 

(3) 


(4) 


The total derivative d( )/dt y x y y, z, and t have the 
standard definitions. The Coriolis parameter / satisfies 
/ = 2ft sin 4> y where ft is the angular velocity of the 
earth and <f> is the latitude. Here p is the density, g is 
the acceleration due to gravity, c p is the specific heat 
for constant pressure, L is the latent heat of conden- 
sation or evaporation, and Q is the condensation or 
evaporation rate. Here Q is a sink for condensation 
and a source for evaporation and sublimation from 
saturated air; otherwise, Q- 0 when q is below satu- 
ration, This is expressed formally by 
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Q = 


dq, 

dt ' 


(5) 


-1 dp 
Ugm fpdy- 


(9) 


Evaporation from the earth’s surface is not explicitly 
included in the foregoing inviscid model. Adding the 
continuity equation to the four preceding equations 
enables the vertical motion to also be computed. For 
large-scale flow the hydrostatic and ideal gas laws com- 
plete the list of equations. 

The transfer of moisture or heat per unit area per 
unit time represents the moisture or heat flux. The 
kinematic flux form is consistent with atmospheric 
measurements of temperature, mixing ratio, and wind 
speed and direction (Stull 1988). For example, the ki- 
nematic moisture flux has components qu , qv y and 
qw y while the same quantities for the kinematic sensible 
heat flux are 8u y 8v y and dw. To find the component 
equations governing the temporal and spatial distri- 
bution of the horizontal kinematic moisture flux, mul- 
tiply ( 1 ) and (2) by the mixing ratio q and rearrange 
by application of the chain rule and (4). (A similar 
process involving the potential temperature yields the 
kinematic heat flux equations described in the Appen- 
dix.) The equations for the horizontal components of 
the total moisture flux are 


djqu) 

dl 


~y + qfo + Aqu, 
p ox 


djqv) 

dt 


~ dfu + Aqv. 

p dy 


( 6 ) 

(7) 


Along with the total derivative, note that each nonlinear 
flux component equation contains a pressure gradient 
term, a term containing the Coriolis parameter, and a 
term that includes the factor A = Q/q . The A 
coefficient 1 is negative for condensation, making Aqu 
and Aqv Rayleigh damping terms. Thus, regardless of 
the sign of qu and qv , the damping term always reduces 
the magnitude of the flux. In contrast, evaporation in- 
duces a forcing response that enhances the magnitude 
of the flux. In tropical regions the term containing the 
Coriolis parameter / is small and may be ignored 
without introducing large errors. Then, changes to the 
kinematic flux depend on the residual of the pressure 
gradient and damping terms. 

Away from the tropics the pressure gradient can be 
expressed in terms of the large-scale geostrophic wind 


1 dp 

Vg **TpTx' 


( 8 ) 


1 Both the moisture and heat flux have a term containing the con- 
densation or evaporation rate. Rayleigh damping reduces the moisture 
flux, while the heat flux is enhanced by diabatic forcing during pre- 
cipitating events. In magnitude, A is significantly larger than the B 
coefficient associated with the forcing in the sensible heat flux equa- 
tions given in the Appendix. 


It should be noted that diabatic contributions to the 
pressure field are implicitly included in the definition 
of the geostrophic wind. In addition, the Rayleigh 
damping term is important, while other indirect and 
feedback diabatic influences on the large-scale flow oc- 
cur, but are not considered in this study. Using the 
geostrophic approximation simplifies (6) and (7), 
yielding 

d( qu) 

~Y t J = -qf(v g -v) + Aqu , (10) 

d(qv) 

—^- = df(Ug-u) + Aqv. ( 11 ) 

The moisture flux is conservative in ( 10) and ( 1 1 ) for 
adiabatic and frictionless conditions, and for geo- 
strophic flow. If the wind components were to remain 
geostrophic during a condensation event, then signif- 
icant reductions in the kinematic moisture flux would 
occur because of the cumulative effects of the Rayleigh 
damping. The exact nature of any optimal local re- 
sponse in the moisture flux, heat flux, energy, or other 
basic atmospheric state is unknown. Scale analysis 
(given below), however, suggests that the large-scale 
moisture flux is nearly conserved. Even though the 
surface evaporation process slowly builds up the at- 
mospheric water vapor, condensation can quickly de- 
plete local supplies unless there exists a mesoscale cir- 
culation to help replenish the local reservoir. Any hy- 
potheses that say there is a tendency toward 
conservation of moisture flux must also include some 
mechanism for generating vertical circulations. 

One way to make the left-hand sides of ( 10) and 
(11) small in magnitude is to alter the wind compo- 
nents u and v so that the right-hand sides are closer to 
being balanced but not necessarily identically balanced. 
It is our hypothesis that the reduction in the large-scale 
horizontal wind components observed during major 
long-lived precipitating events help preserve the large- 
scale moisture flux. This reduction in velocity speed 
decreases the change in the synoptic-scale kinematic 
moisture flux that would otherwise occur if the wind 
speed remained unchanged. As will be shown in the 
example and applications in sections 3 and 4, this 
modification also enhances the longevity of the dis- 
turbance and minimizes loss in the moisture flux by 
creating a mesoscale vertical circulation that helps re- 
store or pump additional moisture into the storm. 
Much of the discussion that follows examines the na- 
ture of the changes in the large-scale horizontal wind 
components that are needed to support our hypothesis. 

A scale analysis of ( 1 ) and (2) shows that for large- 
scale conditions the pressure gradient and Coriolis force 
terms are larger in magnitude than the contribution 
from the total derivative (Holton 1979). In a similar 


1995 


MONTHLY WEATHER REVIEW 


Volume 121 


Table I. A scale analysis* of the terms in the 
u moisture wind [Eq. (8)]. 



diou) 

dt 

qfi> 

qdp 

p dx 

Aqu 

Scaling 


qltlU 

q A p 

U'q 

Magnitude 

(extratropics) 

L 

q X 10"* 

q X 10‘ 3 

P L 

q X 10" 3 

U 

<q X 10~ 3 

(tropics) 

q X i0- 4 

X 10" 4 

q X 1(T 3 

<q X 10~ J 


* U ~~ 10 m s‘\ A ~ \ L ~~ 10 6 m, L c (precipitation scale) 

10 5 m, Ap — 10 3 Nm“ 2 . 


analysis for (6) and (7) the term containing diabatic 
processes is an order of magnitude more important 
than the contribution from the total derivative, pro- 
vided A is scaled with Q proportional to 10“ 7 g g -1 s" 1 , 
a magnitude common to mesoscale precipitation 
events. Table 1 illustrates our findings for large-scale 
circulations containing a long-lived mesoscale diabatic 
process, of length scale L c , for tropical and nontropical 
situations. In Table 1 the derivative term yields a mag- 
nitude U 2 q/L , while the diabatic term gives U 2 q/L c , 
where Q is scaled by Uq/ L c and L c - 10 5 m. In tropical 
regions the Coriolis parameter is reduced; thus, the 
pressure gradient and the damping terms dominate 
provided A is sufficiently large. This scale analysis sug- 
gests that the large-scale fluid flow is nearly balanced 
by the right-hand sides in ( 6 ) and ( 7 ) . This is consistent 
with our stated hypothesis. Clearly, different scaling 
will provide different magnitudes; for example, short 
duration or microscale diabatic processes are of much 
less importance to the large-scale flow. However, the 
Rayleigh damping term remains important provided 
LJL '—0.1. Convective parameterizations assume a 
similar ratio between grid and diabatic scales. However, 
the total derivative in ( 10) and (11) becomes important 
as the scale size is reduced. 

In this study ( 10) and ( 1 1 ) or (6) and (7) are eval- 
uated diagnostically. The right-hand sides of these 
equations are assumed to be nearly in balance, but 
some contribution from the total derivative is allowed. 
A simplified analysis is used to trace how the noncon- 
servative total derivative contributes. As an example, 
assume that the deviation from conservation is known 
to be a small fraction of the flux; then for nonzero 
values of q y u, and v there exist coefficients a u and a v 
such that 


d{qu) 

dt 


a u qu. 


d{qv) 

dt 


a v qv . 


( 12 ) 

(13) 


This is a generalization of the approximation for the 


partial time derivative used in Gill ( 1980) and Matsuno 
( 1966). To get a grasp of the magnitude of the coef- 
ficient assume that the flux qu experiences a complete 
reversal of sign for the synoptic flow within a time pe- 
riod of A/; then the coefficient a u would have magni- 
tude 2/A/. Generally the magnitude is less than this. 
The coefficients a u and a v , however, are locally func- 
tions of time and the independent spatial variables. 
For long-lived large-scale diabatic situations these coef- 
ficients vary slowly in space. Within mesoscale distur- 
bances they should vary in a pattern consistent with 
the disturbance. For example, the convection may 
moisten the upper atmosphere, giving a u and a v a pos- 
itive sign, while drying lower levels making the coef- 
ficients negative. Equations (12) and (13) allow a cli- 
matological estimate of the nonconservative property 
of the horizontal kinematic moisture fluxes to be easily 
determined from model calculations. A more complete 
solution procedure that makes no simplifying assump- 
tions will be presented elsewhere in an independent 
study. 

Substituting (12) and ( 1 3 ) into ( 1 0 ) and (11) gives 


-jv g +fo + u(A - a u q) = 0, ( 14) 

fu g - fu A- v(A - a v q) = 0. (15) 


Note that A , a u q y and a vQ couple the two algebraic 
equations together. Dividing (14) and ( 15) by / and 
solving the coupled system for u and v yields 


- it M v v g 
I + M U M V 9 


(16) 


v g M u Ug 
1 + M U M V ’ 


(17) 


Here the subscript m denotes the moist solution. Note 
that u m - u g and v m = v g when M u = (A - a u )f~ l and 
M v * (A — at v )f~ ] are set to zero. Here A is pegative 
for condensation, so a positive value of a u and a v would 
enhance the magnitude of M u and M v . For some cir- 
cumstances — for example, for large hourly condensa- 
tion rates in a hurricane — the magnitude of the non- 
dimensional number R, * A f~ ] can exceed unity. A 
condensation rate of Q ^ -0.8 X 10~ 7 g g _l s _l yields 
Ri = -0.408, provided / = 10~ 4 and an average q 
satisfies q - 2 g kg _l . Magnitudes for R x of 0.2 and 
less, however, are probably more common for large 
areas. 

For tropical flow the moist wind is computed using 
(6), (7), ( 12), and (13). This yields 



(18) 

(19) 


provided neither a u - A nor a v - A is equal to zero. 
The perturbation in the horizontal velocity is found 
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using the pressure gradient computed from the diabatic 
heating. If the total derivative term; that is, in (6) and 
(7) is small, then the perturbation flow induced by the 
diabatic process is antitriptic; that is, there is a near 
balance between the Rayleigh damping term and the 
pressure gradient term. An example of antitriptic flow 
is found in Raymond ( 1986), illustrating topograph- 
ically induced mesoscale circulations in barotropic flow 
around an isolated mesoscale obstacle. Similar circu- 
lations containing positive and negative couplet pairs 
of convergence-divergence and cyclonic-anticyclonic 
flows are produced by this antitriptic mechanism. 

For exlratropical flow an examination of ( 16) and 
(17) expanded in powers of Rj provides a clearer pic- 
ture. Assuming that u m — u L + and v m = v L 
+ where u L and v L are the horizontal wind com- 
ponents in the large-scale flow, and assigning all terms 
involving Q to the perturbation, yields for the pertur- 
bation and large-scale flow the expressions 


Mpcrt = RlVg - R]u g + 


u L = H, 


2., , ^i(«u + a v )u g 


f 


+ • • •, (20a) 


f P 


+ •••, (20b) 


Vpcn = ~R>Ug - *?D S + * ,(gu + 


/ 


+ 


(21a) 


a u u g a u a v v g 

V L ~ V g + — + • • • . 


(21b) 


Note that the large-scale flow is nearly balanced if 
\a v \/f < 1 and | a u \!f< 1. Also the series in (20b) 
and (21b) converge rapidly, provided these ratios are 
much less than 1 . When this is not the case, this series 
expansion should not be used. 

In the Northern Hemisphere, given Rj < 0 (con- 
densation) with u g > 0 and v g > 0, (20a) and (21a) 
imply that the magnitude of a westerly flow is reduced, 
since u g /u m > 1, while the magnitude of a southerly 
flow is enhanced, since v g /v m < 1. Repeating this 
thought experiment with an easterly flow leads to the 
reverse conclusions. In the Southern Hemisphere, u g 
> 0 implies that a northerly flow is enhanced because 
the coriolis parameter changes sign making R, > 0. As 
will be shown in the hurricane example given below, 
the diabatic perturbation terms [(20a) and (21a)] add 
significantly to the local convergence. 

The diabatic contribution alters the speed V m 
= (wj,l fm) l/2 - Writing the definition of speed in 
terms of the large-scale wind and the diabatic pertur- 
bation gives 

V m = [l'l+ 2u L U pen + 2Dz.Dp.rt 

+ (Wpcn) 2 + (Dpert) 2 ] l/2 . (22a) 

Here V L = (u[ + d!) i/2 . Approximating the large- 


scale flow by factoring out V L but replacing all occur- 
rences of V i in the denominator by V g yields 



(22b) 


where V g = (wf + d|) i/2 . Substituting from (20) and 
(21 ) into ( 22b) gives 


V — V 

’ m ' g 


\ - R 2 ,+ 


2 , 2R| (Ot U 4* tty) 


/ 


__ 2R|(tt^t/g 4~ cx v v g ') 2R\U g v g (ot v tt w ) 


fV\ fV\ 


2R)U g v g (a u + ofy)(tt 0 — tt w ) 

Tv\ 


11/2 


+ 


(22c) 


Note that V m = V g { 1 - R?) 1/2 when a u = tt u = 0. All 
terms in the brackets except for the first two tend to 
nullify each other or make small contributions pro- 
vided |tt v | //< land | tt w |//< 1. Thus, (22c) implies 
that the speed associated with moist cyclonic activity 
is less than the geostrophic value. 

For comparison purposes, note that the relationship 
between the geostrophic and gradient wind speed V# 
is 



where R is the radius of curvature (Holton 1979). For 
cyclonic situations (/R > 0) the ratio V v /V g is less 
than 1. Thus, the speed is reduced from geostrophic 
values by the gradient wind primarily because curvature 
effects are included. In regions having condensation, 
the rate needed to sustain a specific curvature can be 
determined approximately by equating (22c) and (23), 
and solving for Q. Approximating the solution by the 
first two terms in (22c) gives 



As an example, if /= 10~ 4 , V/fR - 0.2, and q = 2 
gkg“ ! ,then Q = -1.1055 X 10~ 7 gg -1 s" 1 . Note that 
Q is proportional to fq . 


3. An example 

Idealized flow in a hurricane is used to illustrate the 
diabatic influence on the wind field. An idealized axi- 
symmetric vortex flow was selected because of its in- 
teresting theoretical properties and because of the ini- 
tialization difficulties encountered when inserting a 
diabatic vortex flow in a numerical model. Following 
the procedure presented in section 2, but for cylindrical 
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coordinates and for axisymmetric flow, allows the moist 
radial wind component u m to be written in terms of 
the tangential velocity v according to 

u m = vM v [\ +Jj (25) 

It is assumed that the contribution from a v f~ l in M v 
is small. Near the vortex core v/fr > 1; thus, u m sim- 
plifies to 

u m =frM Vi (26) 

while far from the core v/fr 1, so 

u m = vM v . (27) 

The divergence is 



In the absence of condensation or evaporation ( Q 
= 0) our solutions satisfy u - 0 unless a v # 0. Thus, 
radial velocities associated with condensation bands in 
hurricanes should be localized and exist in addition to 
a dynamically induced large-scale convergence asso- 
ciated with a v ¥= 0. These bands in the radial velocity 
are seen in both numerical model simulations and in 
observations of hurricanes (Willoughby 1979). 

With a balance between the centrifugal and pressure 
gradient forces that characterizes intense atmospheric 
vortices, cyclostrophic flow leads to a tangential velocity 
defined by 



Therefore, solving for the moist tangential velocity v m 
= u, in the cylindrical moisture flux equations, using 
a procedure analogous to that developed in earlier sec- 
tions, gives 

v„ + frv m +fru m M u = v 2 . (31) 


For v/fr > 1, the divergence [(28)] simplifies to 

^ + ^ = 2/A/„ + ^, (29a) 

dr r dr 

while when v/fr 1 the expression is 

£ 5 + « ! = *.(i’ + SU£^. (29b) 

dr r \dr rj dr 


In (29a) the sign depends upon M v and its radial vari- 
ation, In (29b) the divergence depends upon the vor- 
ticity (dv/dr + v/r) and the radial variation in M v . In 
both (29a) and (29b) M v < 0 contributes toward the 
convergence, provided the magnitude of M v is increas- 
ing radially outward and v > 0. Positive vorticity, av- 
eraged over the hurricane inner core, of 10" 2 $ _1 is 
common. Thus, as indicated by (29), substantial con- 
vergence is generated by diabatic processes when M v 
is sufficiently large in magnitude and of negative sign. 
Negative vorticity with M v < 0 and decreasing outward 
implies divergence. Evaporation near the cloud tops 
also encourages divergence. 

To express the radial velocity as a function of dis- 
tance from the vortex center, it is necessary to know 
M Vi or Q when a v f~ ] is small, and its radial distri- 
bution. Near the vortex center a constant value for M v 
would imply, following (27), a flow that is proportional 
to r (solid rotation), while if M v were linear in r, then 
u would be proportional to r 2 . Far from the vortex 
center the radial component [(28)] is proportional to 
the tangential component v and to M v . In the outer 
reaches of a hurricane the radial distribution of v de- 
creases approximately as r“ ,/2 (Riehl 1963). Thus, the 
radial variation in u m is considerable but depends a 
great deal upon the radial variation of the net conden- 
sation rate Q . 


Ignoring the dependence of u m on v m yields a moisture- 
modified gradient wind expression 


2 ~ 2 


j _ 4a m A/ K + 


Av] 


fr (fr ) 2 


1/2 


(32) 


Retaining the dependence of v m on u m developed above 
[(25)] requires solving the cubic equation 

Vi + 2 frvi + [( 1 + M v Mu(fr) 2 ~ v 2 ]v m = frv) 

(33) 


for tv Note that (33) remains cubic even when M v 

- M u = 0, implying than an extraneous real root 2 has 
been introduced by the algebra. 

Results are presented in Fig. I to illustrate the dif- 
ferences between (32) and (33). The axisymmetric 
tangential flow is taken from the analytical expression 
given in(6)ofFiorino and Elsberry ( 1989), with max- 
imum velocity K mtx = 35 m s“‘ . In their formula for 
the tangential wind the radius of maximum tangential 
winds is set to r max = 100 km and the parameter b 

- 0.96. For simplicity, in our experiments this radial 
profile is assumed to be in cyclostrophic balance, and 
is illustrated by the thick line in Fig. 1. Adding the 
influence of the Coriolis parameter (gradient wind) for 
a dry case [u m = 0, (32)] produces reductions in the 
tangential velocity of approximately 4 m s _I at r max , 
with larger values outside the radius of maximum 
winds. Willoughby ( 1990) indicates that the gradient 
balance condition closely approximations the axisym- 
metric tangential winds observed in hurricanes. How- 
ever, it should be remembered that measured pressure 


2 For the v c used, ( 33 ) has three real roots, two of which are negative. 
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Fta.-l. Idealized axisymmetric hurricane cyclostrophic gradient 
and moist tangential veIocities(m s-')asa function of radius. Speed 
and radial component also illustrated for the moist wind. 


gradients already contain diabatic contributions. This 
could explain the discrepancy between these observa- 
tional findings (Willoughby 1990) and the theoretical 
and numerical findings that intense vortices react 
specified heat sources (Hank and Schubert 
1986; Delden 1989). Note that our moist wind solu- 
tion, as illustrated by the divergence in (29), is also 
dependent upon the intensity of the vortex. 

. . To (33) for v m , the nondimensional values of 
M v and M u (Fig. 2) are assigned to be proportional to 
v . c . so ^ at each is -2 at r max . This large magnitude for 
and M u is used to dramatize the response of the 
drag terms. The roots of (33) give one real positive 
solution for y m that can be used to evaluate (25), yield- 
ing a radial inflow of approximately 12 ms" 1 . As ex- 
pected the diabatic influence reduces the speed and 

inn S t ™ US °f. 7 iaxim “m tangential velocity from 
tUO to 80 km, while generating the maximum mag- 
nitude of the radial velocity at r = 140 km. The max- 
imum reduction in the tangential velocity component 
is about 8ms .but the overall reduction in the speed 
is less since a strong radial inflow is present. Considering 
that our estimates of M v and M u are inflated in mag- 


nitude and in areal coverage, we can probably infer 
that the tangential velocities are not greatly altered by 
the diabatic drag term. 

Magnitudes of M v and M u dictate the radius of max- 
imum heating (RMH). Figure 2 shows a curve for M 
= M V = M u where the RMH is at 100 km, and another 
M that is positive between 0 and 50 km, illustrating 
possible modifications (exaggerated) because of inner 
core evaporation. These M's are used in the calculation 
of u m to show the importance of the divergence as a 
function of the RMH, as illustrated in Fig. 3. In every 
case the balanced solutions show that the greatest con- 
vergence is located interior to the circle with the pre- 
scribed RMH. Whether this would result in a change 
in the location of the RMH depends on the availability 
of moisture and other factors. However, in these cal- 
culations the radius of the maximum cyclostrophic 
wind was held fixed at 100 km. Otherwise, little change 
occurs in the magnitude or location of the maximum 
value of v„ as a function of RMH. 

Observations or results from model simulations of 
intense vortices are needed to complete the evaluation 
of the moisture-modified horizontal winds and the 
evaluation of the diabatic and total derivative terms. 
The contribution from the total derivative would be 
expected to vary spatially, but the pattern should be 
consistent between similar types of situations and 
should change in time somewhat slowly. 

4. Applications 

Cumulus initialization and parameterization of long- 
hved convection in the extratropics are areas that could 
benefit from the knowledge contained in ( 16) (17) 
(18), and (19). The following discussion highlights 
some simple tests to investigate this potential in cu- 
mulus parameterization. These tests are presented to 
illustrate possible applications and do not represent a 
completely definitive study of the topic. 

Cumulus parameterization schemes— for example, 
Kuo (1965, 1974) — estimate subgrid-scale moist con- 



Rc. 2. Radial profiles of M illustrating (exaggerated) radial 
distribution of the condensation or evaporation rate. 
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Fig. 3. Radial distribution of the moist divergence (s' 1 ) 
as a function of M . 
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vective activity using the instantaneous state of the grid- 
scale variables and moisture or cloud properties. This 
information is then used to modify the large-scale tem- 
perature and moisture fields. Because cumulus clouds 
contain a large vertical component of motion, relative 
to their environment, they also are responsible for sig- 
nificant vertical momentum transport. Several schemes 
based on cumulus friction concepts have been intro- 
duced to adjust the grid-scale momentum (Ooyama 
1971; Shapiro and Stevens 1973; Schneider and Lind- 
zen 1967; Shapiro and Stevens 1980; Zhang and Cho 
1991). These approaches introduce a cloud velocity 
that-requires an additional determination, usually from 
a cloud model. In contrast, the terms in (20a) and 
(2 la) that represent approximate changes in the large- 
scale wind introduce no cloud-scale dependent vari- 
ables but require calculation of the geostrophic wind. 

To test the momentum adjustment suggested by the 
perturbation terms in ( 20a) and (21a), they were added 
to the Kuo cumulus parameterization scheme in the 
local [Cooperative Institute for Meteorological Satellite 
Studies (CIMSS)] version of the Australian Bureau of 
Meteorology Research Centre’s (BMRC) finite-area 
model (Leslie et al. 1985). Thus, Wnew W 0 id "h Wpcrt 


where (20a) defines similarly, v ncw = u <*, + 
where (21a) defines Our configuration of this 
semi-implicit regional model has 15 sigma levels. A 
12-h control forecast, 0000-1200 UTC 22 April 1988, 
with 150-km horizontal resolution and 600-s time steps 
was performed using the standard Kuo cumulus pa- 
rameterization scheme ( 1965) without any momentum 
adjustment. In the test forecast the first two terms in 
(20a) and (21a) are used to modify the horizontal wind 
components u and u, with the added requirement that 
| R\ | < 0.5. (Smaller magnitudes for R\ may be more 
acceptable.) This large magnitude for R, was used to 
highlight the differences between the two schemes in a 
short forecast. Even though the convection may vary 
slightly in location and intensity between model time 
steps, it is still assumed that the convection is long 
lived. Spatially isolated convection will not result in 
introducing numerical noise because such noise is re- 
moved efficiently by the high-order low-pass tangential 
or sine filters (Raymond 1988; Raymond and Garder 
1991). In Fig. 4 the differences between the forecasts 
for the 500-mb mixing ratios (g kg' 1 ) are superimposed 
upon the infrared satellite image. Note that small in- 
creases in the mixing ratio occur over cloudy regions, 



Fig. 4. Difference (X 10) in the mixing ratio(g kg“ 1 )fora 12-h forecast (0000-1200 UTC 22 April 1988) 
with and without the perturbation wind [Eqs. (20a) and (21a)] application in the cumulus cloud param- 
eterization scheme. Differences are superimposed upon the VAS infrared image. 
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but decreases occur over regions associated with de- 
scending vertical motion. In this figure the band of 
cjpuds to the north and west is associated with a cold 
front, while the band over the Qulf coast is high clouds 
streaming in over Mexico from the Pacific Ocean. 
(These bands are separated by a strong wind jet.) Our 
tests show that the. first term in (20a) and (21a) is 
responsible for most of enhanced vertical circulation. 

The convective rainfall from the cumulus scheme 
with momentum adjustment was slightly less than that 
produced by the traditional Kuo procedure because 
the restriction | R, | < 0.5 limited the amount of con- 
densation. This made the parameterization more re- 
alistic since too much condensation was occurring too 
rapidly in some localities with the standard scheme. It 
is known — for example, Krishnamurti et al. ( 1980) — 
that the original Kuo (1965) scheme gives excessive 
moistening of the atmosphere. Sundqvist et al. ( 1989) 
introduced a special factor to govern the vertical dis- 
tribution of convective heating to help correct that 
problem. Nevertheless, a maximum increase in the 
nonconvective rainfall of 4 mm in 12 h was also ob- 
served. This increase (Fig. 5) occurred because of the 
enhanced vertical circulation that resulted from the 


application of the momentum adjustment [(20a) 
and (21a)]. 

Initialized National Meteorological Center global 
analyses valid over North America between 1 200 UTC 
9 March and 1200 UTC 12 March 1988 were used as 
initial and boundary conditions in a 72-h forecast made 
by the CIMSS regional model using 190.5-km resolu- 
tion on a 41 X 37 horizontal grid with a 720-s time 
step. The control forecast uses no momentum adjust- 
ment. In the test case the momentum adjustment uses 
the first two terms given in (20a) and (21a) with the 
restriction that | R t | < 0.25. Also, the sixth-order im- 
plicit tangent filter (Raymond 1988) is applied with a 
filter factor of either < = 0.05 or t — 0.64. The latter 
value gives moderate to heavy filtering. The u velocity 
component rms differences between verifying initial- 
ized analysis and forecast are shown as a function of 
time in Figs. 6a-c. The momentum adjustment at 
sigma level 3 (approximately 230 mb) results in a small 
improvement in the statistics except at 24 h. Increasing 
the filter factor from 0.05 to 0.64 has a larger impact 
than the adjustment, primarily because all grid points 
are involved. In contrast, the maximum number of 
occurrences for cumulus parameterization is no more 



Fig. 5. Same as Fig. 4 but illustrating the 12-h difference in the total rainfall (mm) 
between the two 12-h forecasts. 
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Time 


Fig. 6. Root-mean-square differences are illustrated for the u ve- 
locity wind component at (a) sigma level 3, (b) sigma level 6, and 
(c) sigma level 10. “Modified” means the momentum adjustment is 
used. These statistics are computed with an epsilon filter factors (ep) 
of 0.05 and 0.64 in the implicit tangent filter. 


than 527 vertical profiles out of about 1 500 possible 
locations (see Fig. 8). At sigma level 6 (Fig. 6b, ap- 
proximately 440 mb) the major improvement occurs 
late in the forecast. Lower levels (approximately 720 
mb) respond best to the momentum adjustment, as 
illustrated in Fig. 6c. In all cases, heavy filtering always 
reduces the contribution from the adjustment. The v 
velocity wind component rms differences are shown in 
Fig. 7. Again, sigma level 10 shows the most improve- 
ment at 60 and 72 h. Upper-level contributions average 
out to near neutral or slightly on the negative side. This 
may reflect shortcomings in the cumulus parameter- 



Time 

Fig. 7. Same as Fig. 6 but for the v velocity wind component 
and with the smaller filter factor. 


ization or in the momentum adjustment process or in 
both. The lack of entrainment or vertical mixing by 
turbulence in the cumulus parameterization is clearly 
a deficiency. How case dependent the adjustment pro- 
cedure is remains unknown, as do estimates for a u 
and a,. 

The number of active cumulus parameterizations 
versus time step is shown for a 36-h period in Fig. 8. 
It shows that cumulus convective activity is slightly 
enhanced early in the forecast and experiences less 2 Al 
and other high-frequency oscillations with the mo- 
mentum adjustment. 



Time Step 


Fig. 8. Number of grid points experiencing cumulus 
parameterization as a function of the time step. 








JULY 1993 


RAYMOND 


2002 


5. Summary 

A diabatically induced term is present in the equa- 
tions describing the moisture flux and heat flux. This 
term has the form of damping (forcing) for conden- 
sation (evaporation). Diabatic processes also influence 
the horizontal flux through the pressure gradient term. 
In this study the influence of the damping term is eval- 
uated, but the characteristics associated with diabatic 
changes to the geostrophic wind are not analyzed. 

This study finds that mild reductions in the large- 
scale wind speed caused by cumulus friction can con- 
tribute to the near conservation of the large-scale hor- 
izontal components of the moisture flux. This is the 
consequence of an establishment of a near balance be- 
tween the pressure gradient, Coriolis, and Rayleigh 
damping terms in the moisture flux equations. This 
balance is identified by scale analysis. This outcome is 
independent of the mechanism responsible for the re- 
duction in wind speed. 

Estimates of diabatically related changes to the hor- 
izontal wind are easily obtained from the balanced re- 
lationship. For realistic condensation rates, these wind 
estimates are consistent with the gradient wind ap- 
proximation. When these estimates are used in a simple 
hurricane model they predict a radial inflow with con- 
vergence, provided the rotation is cyclonic. In contrast, 
divergence is associated with anticyclonic rotation and 
evaporation. When the predicted wind perturbations 
are used as a momentum adjustment in a Kuo cumulus 
parameterization scheme they improve the 72-h fore- 
cast rms u velocity component error statistics and en- 
hance mesoscale vertical circulations. In the tropics 
any movement toward conservation of moisture flux 
enhances an antitriptic relationship, cementing a near 
balance between a Rayleigh damping term and the 
pressure gradient term. Note that it is necessary to have 
the condensation before the wind-related processes can 
be activated. In this way this process is similar to con- 
ditional instability of the second kind (CISK). 

Our equations (20a) and (21a) are approximations 
of the influence from one diabatic source on the wind. 
A more rigorous study using the vorticity and diver- 
gence concept will be described elsewhere. The terms 
in the flux equations are also being analyzed using me- 
soscale model calculations. This is the subject of con- 
tinuing study. 
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APPENDIX 


The Sensible Heat Flux Equations 

The equations describing the temporal and spatial 
distribution of the horizontal components of the sen- 
sible heat flux can be constructed from ( 1 ), (2), and 
(4). They are 


d(0u) 

dt 


6 dp 

+ Bjv — B6u y 

pdx 


d(9v) 

dt 


— ~ — 6fu — B6v. 
pdy 


(Al) 

(A2) 


In these equations the damping coefficient B = LQ / 
( Tc p ); here, T is temperature, c p is the specific heat at 
constant pressure, and L is the latent heat of vaporiza- 
tion. For condensation processes B is negative but 
smaller in magnitude than the coefficient A in (8) 
and (9). 
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